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We describe a variational procedure for calculating the energy of an electron gas in which the 
long-range Coulomb interaction is truncated by the screening effect of a nearby metallic gate. We use 
this procedure to compute the quantum capacitance of the system as a function of electron density 
and spin polarization. The accuracy of the method is verified against published Monte-Carlo data. 
The results compare favorably with a recent experiment. 



The electron gas with 1/r Coulomb interactions is a 
fundamental reference system. It provides both a testing 
ground for many-body theory and a microscopic input 
for the density functional theory (DFT), which is the 
main tool of electron structure calculations pQ. However, 
there are compelling reasons to study electrons with in- 
teractions other than Coulomb. One is a recent interest 
in mixed DFT schemes, where the Coulomb potential 
is split into a long-range part, to be handled by other 
techniques, and a short-range part, to be treated within 
the usual DFT. This allows one to capture long-range 
van der Waals interaction effects but necessitates recom- 
puting the exchange-correlation energies for the trun- 
cated Coulomb potential [2]. Another and more direct 
motivation comes from the physics of low-dimensional 
systems, in which the bare interaction is often modified 
by the environment. One example is a two-dimensional 
(2D) electron gas positioned a small distance d away from 
a metallic gate. The gate creates an image charge for each 
electron, leading to the interaction law 
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where k is the dielectric constant of the medium. At 
r >• d, this potential rapidly decays: V(r) ~ 2e 2 d 2 /nr 3 . 
Therefore, electron correlations at low density n <C 1/d 2 
are very different from those in the absence of the gate. 
Previously, these correlations have been treated using a 
semi-classical theory in which electrons are assumed to 
form a crystal [3l H]. Such an approximation underes- 
timates the energy but it is legitimate in the range of 
densities a 2 /d 4 <C n <C 1/a 2 , where a — h 2 K/me 2 is the 
effective Bohr radius. Outside of this range, a different 
approach is needed. 

Generally speaking, the most reliable results for the 
electron gas have been obtained by quantum Monte Carlo 
(MC) simulations. While some non-Coulomb potentials 
have been examined [5. in applications to 3 He and nuclear 
matter in three dimensions (3D), no simulations for the 
interaction law V(r) in 2D have been reported. Below 



we demonstrate a variational method for calculating the 
total energy E per electron, which is accurate yet simple 
to implement for any truncated Coulomb potential. Our 
approach is to treat the electron charge e in e 2 /r as an 
adjustable parameter and use the corresponding ground 
states as trial states for the system with the desired in- 
teraction law V(r). Since the 1/r interaction does not 
have an intrinsic lengthscale, the variational energies can 
be obtained from the already existing Monte Carlo data 
by a suitable rescaling. 

Experimentally, the regime nd 2 < 1 has been probed in 
several recent studies [6H9] . Our theory can be compared 
with the one |6] done on a LaAlOs/SrTiOs heterostruc- 
ture with d = 4nm. In Fig. [IJa) we present the results 
for the effective capacitor thickness d* , which is related 
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FIG. 1. (Color online) (a) The effective capacitor thickness 
d* in units of d as a function of n 1//2 d for the spin unpolar- 
ized case, with and without the gate screening effect (thick 
black and thin red lines, respectively). The dashed line is the 
previous classical theory [4]. Symbols are the experimental 
data [6]. (b) The ratio of d* for the spin unpolarized (£ = 0) 
and polarized (£ = 1) cases. 
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to E(n) by 

Aire 2 dn Aire 2 dn 2 ^ ^ ^ ' ' ^ ^ 

where p{n) is the electrochemical potential. This quan- 
tity d* can be extracted from the measured differential 
capacitance of the electron gas. Our results are in a good 
agreement with the experiment [5]; however, the gate 
screening effect becomes important only at the lowest 
measured densities [see Fig. [if a)]. Below we first describe 
how these results have been obtained and then show that 
our method captures the essential physics of the problem 
at both low and high densities. Next, we check that our 
solution satisfies the virial theorem and verify the accu- 
racy of our method on a 3D example. Finally, we discuss 
experimental implications. 

Variational procedure. The trial states for our vari- 
ational calculation are the ground-states of the 2D elec- 
tron gas with the e 2 /r interaction on a neutralizing back- 
ground. The total energy per electron E, the interaction 
energy E^ int \ and the pair distribution function (PDF) 
g(r) of this system are known to scale with n and with the 
dimensionless parameter r s = 1/V irna 2 . For example, 

E = E f(r s ), E^=E f^(r s ) (no gate) . (3) 

Here Eq — irh 2 n/m is the Fermi energy of a noninter- 
acting 2D Fermi gas and f(r s ) is a dimensionless func- 
tion pp. By the virial theorem, f^ lnt \r s ) = r s f'(r s ). 
Similarly, we can write g{r) — G{kpr,rs), where G is 
another dimensionless function and kp = is the 

Fermi momentum. We rely on the fact that / and 
G have been computed by MC techniques and fitted 
to analytical expressions 10-12] over the broad range 
1 < r s < 40. This enables us to use e 2 , or equivalently, 
r s , as a variational parameter that labels the trial wave- 
function. We denote this parameter p s to distinguish 
it from the physical r s . Large p s describes a strongly- 
correlated, crystal-like arrangement of electrons [3j 0], 
while p s = corresponds to the free Fermi g as, as 
in the Hartree-Fock (HF) approximation. The optimal 
p s is to be found by minimizing the variational energy 
E v&T (n,p s ) = Eial n \n,p s ) + Eifr\n 1 p s ). We expect 
p s < r s because screening by the gate reduces electron 
repulsion and weakens electrostatic correlations. 

For a given trial p s , the kinetic energy E^ n \n, p s ) is 
obtained by subtracting the interaction energy from the 
total energy for the 1/r interaction: 

E^\n,p s )=E [f(ps)-Psf'(Ps)}. (4) 

By making use of the definition E = E (l + ( 2 )/2 + E x + 
E c , where E x is the exchange energy, E c — E f c (r s ) is 
the correlation energy, and < C < 1 is the spin polar- 
ization, we can rewrite Eq. Q as 

E&=Eo(j^- Ps f'My (5) 
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FIG. 2. (Color online) (a) Optimal p s vs. r s . Solid (dashed) 
lines correspond to £ = (£ = 1). The labels indicate the 
values of d/a. The inset depicts the maximum value of p 3 for 
C = (circles) as a function of d/a (line), (b) Pair distribution 
functions g(r) for the optimal trial states at r s — 100. 

The function f c (p s ) is taken from Ref. [H] The interac- 
tion energy E^ar*' is computed by numerical quadrature 
of the potential V(r) weighted by the PDF, 

oo 

Et?(n,p s ) = \Jv G(x,p s )xdx - £^ . (6) 

o 

(The last term accounts for the interaction of each elec- 
tron with its image.) For G, we use the parametrization 
of Ref. 1121 The minimization of E vaT over p s is done nu- 
merically. Finally, d* is computed according to Eq. 

We performed these calculations for d/a ranging from 
4 to 32 and for two values of the spin polarization, ( = 
(unpolarized) and 1 (fully polarized) [T3"] . The optimal 
p s as a function of the physical r s is shown in Fig. [2^ a) 
for C = 0, 1. At large electron density (small r s ), p s re- 
mains close to r s because the average spacing between 
electrons is much smaller than d, so that the gate plays 
a minor role. As r s increases, p s attains a maximum 
value p s ~ d/a, and then drops off. This is because 
at r s > d/a the interaction between electrons becomes 
effectively short-ranged, so that electrons lose their elec- 
trostatic correlations as the system becomes more dilute. 
For f = 1, the drop of p s at r s > d/a is fairly rapid. 
In contrast, in the unpolarized case, p s exhibits a broad 
plateau before it also collapses at very high r s [far beyond 
the range shown in Fig.[2|a)]. 

In the remainder of this paper we show that our results 
withstand several tests: (i) they correctly reproduce the 
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asymptotic behavior at large and small n, (ii) they obey 
the virial theorem, and (iii) when generalized to 3D, they 
yield a good agreement with MC simulations [2]. 

Low and high density asymptotics. Because of 

its fast 1/r 3 decay, V(r) belongs to the universality 
class of short-range potentials. At low n, the energy 
per particle coincides in the first approximation with 
that of the free Fermi gas, except that it is shifted by 
— e 2 /4«:c?, cf. Eq. ([6]). This is reproduced by our calcula- 
tion since p s (n — >• 0) = 0. This fact further implies that 
d* = d*(C = 0) and d* =d*(( = 1) tend to a/4 and a/2, 
respectively, as n — > [see Fig. [TJa)]. 

The leading interaction correction to the free gas limit 
comes from two-body collisions, which are parametrized 
by the effective hard-core radius [H] b = 2e 2l d 2 /a, where 
7 w 0.5772 is the Euler constant. In the absence of spin 
polarization (£ = 0), the perturbative result for the elec- 
trochemical potential reads [15] u = E [l + 2v + iv 2 (l — 
ln2)], where [16] 1/v = —27 — ln(7m6 2 ). Accordingly, 
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ln(a 2 /nd 4 ) - 6.30 
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This expression is expected to apply at v < 1/2, which 
corresponds to extraordinary low densities, e.g., no? < 
10~ 6 for d = 4a. Indeed, our variational method shows 
that the departure of djj from its zero-density value of 
a/4 develops abruptly as a function of n within a narrow 
interval n x l 2 d < 10~ 3 [see Fig.flja)]. Of course, at such 
n the energy of real electron systems is dominated by 
disorder [T7] , which we do not consider here. 

The abrupt growth of dp caused by the logarithmic 
correction in Eq. ([7| reflects the low-energy behavior of 
the s-wave scattering phase-shift in 2D. It can be con- 
trasted with the gradual increase of d\ for the polarized 
gas, where p-wave scattering dominates. Here the Pauli 
exclusion between like-spin electrons ensures that elec- 
tron pairs do not approach each other too closely and the 
short-range repulsion is more easily satisfied. As a result, 
in the polarized gas, we have d\ — (a/2)[l + 0(nb 2 )]. 

We can also consider the opposite limit, n ^> 1/d 2 , 
where to the leading order p s = r s and so the kinetic 
energy is unaffected by the gate. The correction AE(n, d) 
to the total energy with respect to the reference system 
without the gate is determined by the interaction energy. 
To compute the latter, we rewrite Eq. ^ as 



= g / V(r)[g(r) 
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The difference g(r) — 1 is appreciable only at |r| < 71 -1 / 2 
and gives — 1/n when integrated over all r. Therefore, to 
the order C(l/d) we can set V(r) = e 2 /(nr) — e 2 /{2nd) 
in the integral, leading to AE = (2ire 2 d/ n)n. The cor- 
responding correction to the electrochemical potential is 
A/i = (4ire 2 d/ n)n, which is the classical relation between 
the voltage \ije and the charge density en of a parallel- 
plate capacitor of thickness d. 



The parameter d* , as defined in Eq. Q, represents 
the effective capacitor thickness, with C = e 2 dn/d/j = 
«/(47rd*) being the capacitance per unit area. In gen- 
eral, d* differs from the geometric thickness d, which is 
often called the quantum capacitance (QC) effect [15] , 
Our preceding calculation shows that at high density the 
gate can modify the QC only to order 1/nd 2 . A simple 
expression for this correction can be derived at 77 3> 1 /a 2 , 
where the HF approximation applies. For £ = 0, we get 
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The QC is represented by the last three terms in this 
equation. They account for, respectively, the kinetic en- 
ergy of the Fermi sea, the exchange energy, and the cor- 
rection to the exchange due to screening by the gate. 
Because of the small numerical factor, the last of these 
becomes important only at n 1 / 2 d < 0.25 [see Fig. [lja)]. 

Finally, the regime of strong correlations induced by 
the 1 /r 3 tail of the interaction [3l H] is realized at inter- 
mediate densities, l/fe 2 ^77^;l/d 2 . The energy of such 
a state is largely insensitive to the spin polarization. In- 
deed, our calculation shows that the relative difference 
of dg and d\ is small at such n. On the other hand, this 
difference becomes significant in the low-density limit, 
where the correlations weaken, cf. Fig.[ljb). 

Variational PDF. Our calculation also provides the 
PDF g var (r) of the screened electron gas [see Fig. [2jb)], 
which has a number of merits. Our PDF is strictly pos- 
itive, unlike those of some approximate many-body the- 
ories pp. One can show that this g va .r{r) and the varia- 
tional energy components satisfy the virial theorem: 
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where P = [p(n) — E(n)]n is the pressure and D is the 
space dimension. (The subscript "var" in P and g is 
implicit.) The variational estimates of E and /1, which 



determine the left-hand side of Eq. ( 10 1 , should be reli- 



able. Therefore, at intermediate values of r, which domi- 
nate the value of the integral on the right-hand side, our 
PDF may be a good approximation. At r much larger or 
smaller than the mean electron spacing, it is probably less 
accurate. Thus, we have g(, ar (0)/g V ar(0) = jfti (ps/r s ) 
instead of the exact cusp condition [T] g'(Q)/g(0) = j^zj- 
Test on a 3D model. The accuracy of our method is 
best verified by comparison with MC simulations; how- 
ever, they are not available for our 2D problem. Instead 
we did such a test for a 3D electron gas with the inter- 
action potential V(r) — (e 2 /r) erfc(/j c r) [5] (relevant for 
mixed DFT schemes). The necessary correlation ener- 
gies and PDF for the standard 3D gas were taken from 
Ref. fT9l The results for the cutoff parameter /i c = 0.5/a 
are shown in Fig. [3] The largest difference between our 
variational estimate and the MC results [5] for the to- 
tal energy per particle is ~ 2%, which is a significant 



4 



S3 



1.06 
1.04 
1.02 




10 



FIG. 3. (Color online) The ratio of the total variational energy 
per particle E and the MC energy -Ezecca [5] as a function 
of r s for a 3D electron gas with the interaction law V(r) = 
(e 2 /r) erfc(/i c r) and fi c = 0.5/a (thick black line). The thin 
red line shows the HF result. 



improvement over the HF approximation. This exam- 
ple also illustrates the capability of our method to treat 
other dimensions and/or interaction laws 20 . 

Discussion. Our results for d* can be verified experi- 
mentally by measurements of the differential capacitance 
between an electron gas and a metal gate. Previously, 
structures with a distant gate, nd 2 > 1, have been stud- 
ied. At low densities, d* was seen to be slightly smaller 
than d [TT1 |2"TH2"5] . In fact, this negative "QC" correction 
arises largely from correlations of classical nature [H • 

Our theory enables us to compute the capacitance of 
gated nanostructures in which quantum and correlation 
effects are not mere corrections. Our results agree well 
with the data of Ref. El Rigorous testing would require 
additional experiments at still lower densities and/or a 
strong in-plane magnetic field (to check the predicted 
spin dependence). Additionally, the energy per particle 
we compute here is relevant for understanding impurity 
screening and therefore transport properties of correlated 
2D electron liquids near metallic gates [25H29] . Finally, 
an intriguing direction for future work is to devise a vari- 
ational method suitable for ultrathin gated structures of 
graphene [7HH1 EM G2] , which is an electron system with 
a nonparabolic spectrum. 
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